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The suppression of large transverse momentum (Pt) hadrons in heavy-ion (A-A) collisions as 
compared to their scaled expectation from proton-proton (p-p) collisions due to the interaction 
of hard partons with the hot and dense QCD medium in A-A collisions is experimentally a well 
established phenomenon. Focusing on leading hadrons produced in hard processes, the medium 
effect appears as energy loss. Beyond that, the question is how the lost energy is redistributed in 
the medium. With increased experimental statistics and most importantly the kinematic range of 
the LHC, studying the properties of full jets rather than leading hadrons is becoming feasible. On 
the theory side, analytic models and Monte-Carlo (MC) codes for in-medium shower evolution are 
being developed to describe jets in the medium. In this paper, expectations for medium-modified 
jet observables, the jet shapes, the thrust distribution and the n-jet fraction, are computed with the 
MC code YaJEM for various scenarios of the parton-medium interaction which all are consistent 
with high Pt hadron suppression data. The computation is done at 20 and 100 GeV jet energy, 
corresponding to probing typical RHIC and LHC kinematics, and the possibility to make an unbiased 
measurement of the observables is discussed. 
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I. INTRODUCTION 

Jet quenching, i.e. the interaction of hard partons cre- 
ated in the first moments of a heavy ion coUision with the 
surrounding medium of strongly interacting matter has 
long been regarded a promising tool to study properties 
of that medium [l|, 0, H, 0, H, Q . The guiding principle is 
to study the changes to a hard process well-known in p-p 
collisions which are induced by the medium. A number 
of observables which essentially probe the modification of 
the leading hadron in a jet are available for this purpose, 
among them suppression of single inclusive hard hadron 
spectra in terms of the nuclear suppression factor Raa 
Q , the suppression of back- to-back correlations [1, [§| and 
also single hadron suppression as a function of the emis- 
sion angle with the reaction plane 10]. Most recently 
also preliminary measurements of fully reconstructed jets 
have become available [ll|. 

Single hadron observables and back-to-back correlations 
are well described in detailed model calculations us- 
ing the concept of energy loss [l3, El, [l3|> i-^- under 
the assumption that the process can be described by a 
medium-induced shift of the leading parton energy by 
an amount AS, followed by a fragmentation process 
using vacuum fragmentation of a parton with the re- 
duced energy. However, there are also calculations for 
these observables in which the evolution of the whole in- 
medium parton shower is followed in an analytic way 
[S [li, [l3- Recently, also Monte Carlo (MC) codes 
for in-medium shower evolution have become available 

[B [n [m, m which are based on MC shower sim- 



ulations developed for e^e 



or p-p collisions, such as 
24| . in which the shower 



'Electronic address: ' trenk®phys.jyuTfi] 



PYTHIA [23| or HERWIG 
evolves in vacuum. In a medium-modified shower evo- 
lution, in contrast to energy loss calculations, energy is 
not simply lost but redistributed in a characteristic way. 

Currently the precise nature of the jet-medium interac- 
tion is not known. In pH we have studied three different 
scenarios including medium-induced radiation and a drag 
force in their effect on various observables in the RHIC 
kinematic range using the im-medium shower MC code 
YaJEM (Yet another Jet Energy- loss Model). The aim 
of the present paper is to extend this study to the LHC 
range and to discuss other observables which probe the 
medium-induced redistribution of energy and are hence 
sensitive to the nature of the parton-medium interaction. 

Observables such as jet shapes allow to study the energy 
redistribution caused by the medium e.g. as a modifica- 
tion of the angular flow of energy inside the jet. Thus, 
in this paper we mainly focus on jet shapes, i.e. the 
transverse flow profile of the jet energy, and their in- 
medium modification. Analytical results for this observ- 
able have already been obtained previously [25, 26]. In 
the medium, there is always an ambiguity if a given par- 
ton is considered to be part of the jet or part of the 
medium. Momentum and angular cuts may provide a 
separation criterion, however in turn introduce a bias 
which affects any given jet observable. Thus, we also 
study the effect of such cuts on the jet shape. In addi- 
tion, we compute other jet observables such as the thrust 
distribution or the n— jet fraction which have been used 
in e+e~ collisions for precision tests of QCD [23]. While 
these quantities may not be easy to observe in heavy-ion 
collisions due to the background medium and the pres- 
ence of multiple hard processes in the same event, the 
comparison with the well established vacuum results and 
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the in-medium results obtained with the different MC 
code JEWEL 18] is nevertheless instructive. 



medium with specified transport coefficient q(r]s,r, 

by 



II. THE MODEL 

The MC code YaJEM is described in great detail in [l^, 
[2lj |. therefore we will only give a brief outline here. The 
aim of the code is to compute the evolution of a QCD 
parton shower following a hard process in the medium. 
In the absence of a medium, the evolution of the shower 
is computed using the PYSHOW routine [28| which is 
part of the PYTHIA package 

In the medium, the main assumption of YaJEM is that 
the kinematics or the branching probability of a parton 
are altered during its propagation. For this purpose, it 
is necessary to make a link between the shower evolution 
which is computed in momentum space and the evolu- 
tion of the medium which is computed in position space. 
We assume that the average formation time of a shower 
parton with virtuality Q is developed on the timescale 
l/Q, i.e. the average lifetime of a virtual parton b with 
virtuality Qb and energy Et coming from a parent par- 
ton with virtuality Qa is in the rest frame of the original 
hard collision (the local rest frame of the medium may 
be different by a flow boost as the medium may not be 
static) given by 
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Going beyond the ansatz of [1^, [20| where we used this 
average formation time for all partons we assume as in 
pli] that the actual formation time Tf, can be obtained 
from a probability distribution 
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which we sample to determine the actual formation time 
of the fluctuation in each branching. This establishes 
the temporal structure of the shower. With regard to 
the spatial structure in terms of spacetime rapidity 77^, 
transverse radius r and angle (f) , we make the simplify- 
ing assumption that all partons probe the medium along 
the eikonal trajectory of the shower initiating parton, i.e. 
we neglect the small difference of the velocity of massive 
partons to the speed of light and possible (equally small) 
changes of medium properties within the spread of the 
shower partons transverse to the shower axis. 
YaJEM provides three different scenarios for the effect 
of the medium on the evolving shower. In the RAD sce- 
nario, the idea is that the virtuality of a shower parton 
can grow during its propagation through the medium. 
Such increased virtuality subsequently leads to additional 
medium-induced branching. In practice, we increase the 
virtuality of a shower parton a propagating through a 



(3) 



where the time Tq is given by Eq. ([2]), the time is 
known in the simulation as the endpoint of the previous 
branching process and the integration dC is along the 
eikonal trajectory of the shower-initiating parton. If the 
parton is a gluon, the virtuality transfer from the medium 
is increased by the ratio of their Casimir color factors, 
3/1 = 2.25. 

In the DRAG scenario which is motivated by results from 
modelling QCD-like A'' = 4 super Yang-Mills theories 
via the AdS/CFT conjecture [23|, we assume that the 
medium exerts a drag force on each propagating par- 
ton. The medium is thus characterized by a drag coef- 
ficient £'(?7<5, r, (/), r) which describes the energy loss per 
unit pathlength. In the simulation, the energy (and mo- 
mentum) of a parton a are reduced by 
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dCDiC). 
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The third scenario has been suggested in [3 and 
is included here for comparison. In the following, it is 
referred to as FMED. Here, the modification does not 
concern the parton kinematics, but rather the branching 
probabilities. In this scenario, the singular part of the 
branching kernel in the medium is enhanced by a factor 
(1 + fmed) where fmed represents an average measure of 
the medium effect for a given path through the medium. 
We assume that all the three relevant parameters 
q, D, fmed can be linked with the medium energy den- 
sity up to a (possibly dimensionful) constant factor as 



q, 



D ~ A:[e(C)]^/'' (cosh p(C) -sinhp(C) cos V). (5) 



and 



K / dC[e(C)]^/''(coshp(C)-sinhp(C)cosV'). (6) 



fmed 



Here, p is the transverse flow of the medium at position 
C and '0 is the angle between flow and the direction of 
parton propagation. For most possible paths of partons 
occuring in a realistic medium, the results in p^.[2l| have 
shown that the precise dependence of e(C) has a small ef- 
fect on the shower evolution. In essence paths through 
an expanding medium can be considered equivalent if the 
line integral along the eikonal path of the shower initiat- 
ing parton results in the same value of AQ^ i^E). We 
will use this scaling law extensively in the following (note 
that it is by definition fulfilled in the FMED scenario). 
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The model assumptions presented here with regard to the 
medium effect on the shower are (for the moment) rather 
quaUtative. Complementary pQCD approaches to the 
medium effect such as the high order opacity expansion 
which can also be solved in a MC scheme [31j are more 
focused on treating the quantum interference between 
subsequent scatterings in the medium correctly and can 
hence be used to test the model assumptions made here, 
or ultimately be envisioned as an extension of the present 
work. 



III. SEPARATION OF JET AND MEDIUM 

For a shower evolving in the medium, there is in gen- 
eral exchange of energy and momentum between what 
one considers shower partons and what one considers the 
medium. Note that within the framework of YaJEM, 
in the RAD scenario the shower gains energy from the 
medium by means of the virtuality increase, in the DRAG 
scenario the shower loses energy to the medium whereas 
the shower energy is conserved in the FMED scenario. 
While this appears surprising at first, it is actually rather 
a matter of book-keeping and defining a distinction be- 
tween shower and medium. 

There is no conceptual way to actually distinguish soft 
shower partons from soft medium partons. Within the 
model framework, an artificial distinction is maintained 
since all partons which explicitly appear are created in 
branchings and hence are considered to be part of the 
shower while the medium appears never formulated in 
terms of partons, but only as an effective influence on the 
shower. It is this artificial separation between shower and 
medium which makes for the surprising above results: In 
the RAD scenario, the possibility of soft partons being 
absorbed by the medium and hence transferring energy 
to the medium is not considered, as a result there is no 
energy transfer from shower into the medium and thus 
the shower energy can only increase. In stark contrast, 
in the DRAG scenario it is assumed that the medium 
can absorb energy without medium partons becoming in 
any way correlated with the shower (for example, in elas- 
tic parton-medium interactions, recoiling medium par- 
tons become correlated with the jet direction). As a re- 
sult, there is no possibility to transfer energy from the 
medium into the shower. Finally, in the FMED sce- 
nario the somewhat artificial assumption is made that 
the medium can modify branching probabilities without 
any explicit energy-momentum transfer, hence the energy 
in the shower remains conserved. 

In a more realistic model in which the structure of the 
medium is resolved, one would define a criterion (say a 
momentum scale) based on which partons are removed 
from the shower and become part of the medium. In 
such a model, all three scenarios would lead to a net loss 
of energy from the shower to the medium through the 
appearance of soft partons in the evolution, in addition 
to possible other mechanisms of energy transfer from jet 



to the medium. Unfortunately we do not yet have such a 
model for the medium, since one goal of the hard probes 
program is to establish the structure of the medium in 
terms of its relevant degrees of freedom and their prop- 
erties. However, lacking such a framework which must 
also be able to account fully for the back-reaction of the 
medium to the energy deposition by the jet to be consis- 
tent, we may still ignore the complications of the medium 
by imposing a momentum cutoff and focusing exclusively 
on the hard part of the jet for which the above caveats 
do not apply. 

Based on hydrodynamical models of the medium, a rea- 
sonable cut to separate the perturbative dynamics of in- 
medium shower evolution from the soft, non-perturbative 
dynamics of the medium is PT„,i„ = 2 GeV at RHIC en- 
ergies and PTmin = 4 GeV at LHC energies. Note that 
such a cutoff introduces a substantial bias on the sample 
of jets which enter the computation of a given observable. 
To gauge this effect, we will in the following discuss all 
observables also without the cut imposed. This is done 
with the explicit understanding that these quantities do 
not correspond to observables, as they each neglect an 
important part of the dynamics between jet and medium 
for which there is no reason to be described by pertur- 
bative QCD. On the other hand, quantities where the 
momentum cutoff is imposed should be in principle com- 
parable with experimental results, although in practice 
complications such as the experimental jet finding strat- 
egy or the presence of multiple hard processes in a single 
event need to be taken into account before making an 
actual comparison. 



IV. THRUST AND n-JET FRACTION 

The overall flow of energy in an event can be traced by 
three event shape observables, namely thrust T, thrust 
major Tmaj and thrust minor Tmm- They involve a sum- 
mation over all final state particles in an event. The 
ALEPH collaboration has measured thrust distributions 
in e+e" collisions j27t in the absence of any medium- 
induced final state effect. In this system, PYTHIA is 
able to account for all observables well, thus establishing 
a baseline for the subsequent discussion. A detailed com- 
parison of PYTHIA with the ALEPH data can be found 
in [27]. 

Thrust is defined as a sum over all final state particle 
three momenta Pj 



T 



J2^\P^■^T\ 



(7) 



It measures how well the final state hadrons are aligned 
in an axis defined by the shower originating partons. For 
T = 1 this alignment is perfect, for T = 0.5 the event is 
spherical and no preferred axis can be identified. Thrust 
major is the projection of all particle momenta on this 
axis 
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FIG. 1: (Color online) Distribution of thrust T (dashed), Tmaj (solid) and Tmin (dash-dotted) for a 20 GeV shower initiating 
quark in vacuum and three different scenarios for parton-medium interaction (see text). Shown is the distribution for all final 
state particles (left panel) and for those above 2 GeV (right panel). Note that the DRAG scenario for the given choice of 
medium parameters renders most events below the cutoff, thus a thrust analysis above the imposed cutoff cannot be done. 
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FIG. 2: (Color online) Distribution of thrust T (dashed), Tmaj (solid) and Tmin (dash-dotted) for a 100 GeV shower initiating 
quark in vacuum and three different scenarios for parton-medium interaction (see text). Shown is the distribution for all final 
state particles (left panel) and for those above 4 GeV (right panel). 



l^i \Pi\ 

while thrust minor sums components of momenta which 
are orthogonal to n and 

rp _ Si |P' " "mil fQ\ 
^ min — I I \ ) 

l^i \Pi\ 

where Umi = nr x n. 

While the sum over all final state particles in e+e~ col- 
lisions can be traced back to a single hard process, this 
may not be so in A-A collisions where the number of 
binary collisions is 0(1000). Thus, even with a cut in 
Pr^i^ imposed, the calculated thrust distribution may 
not correspond to the observable ones. 



In Fig. [T] we show the distributions of T, Tmaj and Tmin 
for showers initiated by a 20 GeV quark for a medium 
path with AQ2^j = 15 GeV^, /\Etot = 15 GeV and 
fmed — 1-5. The proportionality between these param- 
eters is chosen based on the criterion that the medium- 
modified fragmentation functions computed for them ap- 
proximately agree in an interval of 0.4 < z < 0.7. This 
is the region of the fragmentation function which is pre- 
dominantly probed when the fragmentation function is 
folded with a pQCD parton spectrum to compute single 
inclusive hadron production. In essence, for this propor- 
tionality between the parameters all scenarios reproduce 
the measured single hadron suppression Raa equally well 
[2l|, and hence the results should in some sense be com- 
parable. The particular choice of ^Qtot — 15 GeV then 
corresponds to a path through the whole RHIC medium, 
i.e. a medium modification close to the maximum which 
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is possible. 

It is evident from the figure that at this kinematics the 
medium tends to make the event more spherical, i.e. the 
distributions widen and move towards 0.5. While the 
effects for all particles are dramatic, the observable dis- 
tribution above the imposed cutoff still shows some mod- 
ification. The effect of both radiative scenarios RAD and 
FMED is qualitatively similar and agrees with what has 
been reported in [ll] using the MC code JEWEL. 
We show the same distribution for a 100 GeV shower 
initiating quark as appropriate for the LHC kinematic 
in Fig. [2l For the sake of comparison, we have chosen 
the medium parameters AQj^^, Ai?tot and fmed identical 
to the RHIC case above. However, since the medium at 
LHC is expected to be somewhat more dense [13], these 
values represent now a fairly typical path through the 
medium rather than the maximal medium effect as in 
the RHIC case. 

It can be noted that without the cutoff the DRAG sce- 
nario acts rather different at higher Pt — here it tends to 
focus the distribution towards T = 1 whereas the radia- 
tive scenarios RAD and FMED still make the event more 
spherical. However, RAD does so only at low Pr whereas 
FMED does so at all Pt — if a cutoff is imposed, it is 
evident that the high Pt part of the thrust distribution 
is not widened in the RAD scenario which is a possible 
distinctive feature between the otherwise rather similar 
RAD and FMED prescriptions. 

Let us now focus on a different observable, the n-jet frac- 
tion. This observable is sensitive to the substructure of 
a jet in terms of clusters of hadrons created by the dy- 
namics of showering and hadronization. It is based on 
the Durham clustering algorithm [33i] . This algorithm 
clusters final state particles based on a distance measure 
between a pair i,j 

y,, = 2mm{Ef, E^){1 ~ cos(%)/ii;2„. (10) 

This yij is a measure for the transverse momentum of a 
softer particle of a pair with respect to the axis defined by 
the harder particle. In each step, the clustering algorith 
replaces the pair with the smallest yij by a cluster with 
energy and momentum given by the sum of the pair's 
energy and momentum. The procedure is repeated and 
particles and clusters are further merged until yij exceeds 
a pre-defined threshold ycut which corresponds to a given 
resolution scale. The number of clusters present at this 
point is called n for this event, given ycut- In averag- 
ing over many events, one finds for each choice of ycut a 
fraction of events with n — 1,2,3,... events — this is 
referred to as the n-jet fraction. 

The same caveats which were relevant for the thrust dis- 
tribution hold for the n-jet fraction: Even with a cutoff 
imposed, one cannot necessarily assume that there is only 
a single hard event in a heavy-ion collision, and the ef- 
fect of multi-jet events must be accounted for before a 
comparison with data can be made. 



We show the 2-jet and 4-jet fraction as a function of the 
resolution scale ycut for a back-to-back pair of 20 GeV 
quarks corresponding to RHIC kinematics in vacuum and 
in medium (with medium parameters chosen as before) 
in Fig. [3l For a coarse resolution (log;^o(ycMt ~ — 1) the 
algorithm picks up two jets corresponding to the back-to- 
back event, when the resolution scale is increased, find- 
ing more clusters becomes increasingly likely. Note that 
some scenarios (especially with the cutoff imposed) sat- 
urate for increased resolution scale. This is caused by 
the relatively small multiplicity in the event, either when 
(as in the DRAG scenario) a large amount of energy is 
transferred to the medium or when a cutoff is imposed 
while soft particle production is enhanced. No increase 
in resolution scale can find six clusters if there are only 
four particles in the final state above the cutoff. 

For RHIC kinematic conditions and with the cutoff im- 
posed, the effect of the medium is to increase the number 
of clusters seen at a given resolution scale. This agrees 
with the previous observation based on the thrust distri- 
bution that the medium tends to make the event more 
spherical — this also implies that the spread in trans- 
verse momenta is increased, and this is chiefly what is 
seen in the plots. 

The 2-jet and 4-jet fraction for LHC kinematics, i.e. a 
back-to-back event of quarks with 100 GeV energy each 
is shown in Fig. [D Here the radiative scenarios RAD 
and FMED show qualitatively the same picture as for 
RHIC — at the same resolution scale they tend to lead to 
more clusters. In contrast to this is the DRAG scenario 
which shows the opposite trend. This agrees with the 
thrust distribution in the DRAG scenario, which show 
that the events become more focused. Combined with 
the fact that overall multiplicity is reduced due to en- 
ergy transfer into the medium, we can thus understand 
the trend of the DRAG scenario. If a cutoff of 4 GeV is 
imposed, the trends qualitatively remain, with the excep- 
tion of the RAD scenario which moves closer to the vac- 
uum result. This again shows that the medium- induced 
changes in the RAD scenario chiefly affect the low Pt 
part of the shower, which was noted before when study- 
ing the thrust distribution with cutoff imposed. These 
findings are qualitatively in agreement with the radiative 
energy loss scenario FMED as implemented in the MC 
code JEWEL [3. 

Unlike other observables such as the nuclear suppression 
factor for single hadrons or the longitudinal momentum 
distribution in the jet studied previously 21], thrust and 
especially n-jet fraction offer thus in principle the oppor- 
tunity to distinguish details of the jet-medium interaction 
even on the level of distinguishing different implementa- 
tions of radiative energy loss if they can be measured 
reliably in A-A collisions. 
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FIG. 3: (Color online) The two (solid) and four (dashed) jet fraction for two back-to-back shower initiating quarks with 20 
GeV energy each in vacuum for three different scenarios for parton- medium interaction (see text). Shown is the distribution 
for all final state particles (left panel) and for those above 2 GeV (right panel). Note that the DRAG scenario for the given 
choice of medium parameters renders most events below the cutoff, thus a clustering analysis above the imposed cutoff cannot 
be done. 




FIG. 4: (Color online) The two (solid) and four (dashed) jet fraction for two back-to-back shower initiating quarks with 100 
GeV energy each in vacuum for three different scenarios for parton-medium interaction (see text). Shown is the distribution 
for all final state particles (left panel) and for those above 4 GeV (right panel). 



V. JET SHAPES 

The jet shape measures the angular distribution of the 
flow of energy transverse to the jet axis. Once a jet axis 
and the particles i associated with the jet have been iden- 
tified experimentally, the integral jet shape given a cone 
radius R is defined as 



Ri 



E,e{R - R,) 



(11) 



where Ei is the transverse energy of particle i, r is an 
opening angle and Ri is the angles of particle i with the 
jet axis. If the jet axis is located at pseudorapidity ry and 



azimuth 0, then Ri = \/ {rji — rj)-^ + 



,R) 



is thus the fraction of energy inside a sub-cone radius r 
if the total energy is contained inside a cone of radius R. 
From the integral jet shape, the differential jet shape can 
be obtained as 



ri^,„t(r, R) 
dr 



(12) 



which is the angular density of energy flow in the jet. 
Since the jet shape depends on the jet energy as jets 
with larger energy are more coUimated, one would like to 
study medium modifications of the jet shape for events 
in vacuum and in medium for the same energy of the 
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FIG. 5: (Color online) Differential energy flux as a function of angle with the jet axis r for jets in vacuum and three different 
in-medium scenarios (see text) with a cutoff imposed for RHIC kinematics (left panel) and LHC kinematics (right panel). 



shower-initiating parton. However, as mentioned before, 
since jet and medium interact, the energy of the shower- 
initiating parton is in general not the energy of the final 
state jet, and moreover one needs some criterion to even 
make a distinction between jet and medium in the soft 
sector. Thus, one needs a handle on the initial energy in 
order to make this comparison. Experimentally, this can 
be done e.g. using 7-jet correlations. 
In order to see that this is indeed a relevant effect and the 
energy of modified jets above the imposed cutoff is not 
roughly the same as the energy of the shower initiator, 
we show the angularly differential flux of energy dE/dr 
for RHIC and LHC kinematics with the relevant medium 
cutoff imposed in Fig. [5l The difference in energy flux 
above the cutoff is more than a factor two between the 
different scenarios even at LHC kinematics. This strongly 
emphasizes the point that it is necessary to get a handle 
on the energy of the shower initiating parton. 
Differential jet shapes for RHIC kinematics, i.e. a 20 GeV 
quark in vacuum and for the different in-medium scenar- 
ios are shown in Fig. [51 for LHC conditions the same is 
shown in Fig. [T] The observed picture is very consistent 
with what was found before: While the medium in the 
DRAG scenario tends to coUimate the jet, the radiative 
scenarios RAD and FMED widen it and lead to energy 
flux at larger angles. However, there is an important dif- 
ference between the RAD and the FMED scenario: In the 
RAD scenario large angle radiation appears only at low 
Pt, hence the observed widening is significantly reduced 
once a cutoff is imposed. This is not so in the FMED 
scenario where the effect persists despite a cutoff. 

VI. JET TRIGGERS AND BIAS 

At least for LHC conditions, one may conclude from 
Figs. [5] and [7] that on average most of the energy flux 
of the jet is within a relatively narrow cone R < 0.5 and 



that even above a cutoff of Pt > 4 GeV the integrated 
flux is substantial, regardless which in-medium scenario is 
studied. Thus, it should definitely be possible to directly 
trigger on medium-modiflcd jets (an energy flux within a 
narrow cone) instead of using 7-jet coincidences. 
The main problem with direct jet identification is to de- 
termine the total energy of the jet and to separate it 
from the background. There are three different contri- 
butions within a given angular region: 1) the (perturba- 
tive) jet itself, 2) medium which is (non-perturbatively) 
correlated with the jet due to the jet-medium interac- 
tion (e.g. through elastic recoil of medium particles, 
Shockwaves, sound mode excitations in a strongly cou- 
pled fiuid medium. . . ) and 3) uncorrelated background 
medium. Since the average energy contained in the un- 
correlated background can be observed event by event, 
and since due to energy conservation the sum of energy 
contained in the jet and the medium correlated with the 
jet must equal the energy of the shower-initiating par- 
ton, one might think of a strategy where a jet is identified 
above a cutoff in Pt and its total energy is determined by 
subtracting the energy content of the uncorrelated back- 
ground and to correct for the jet energy below the cutoff. 
This procedure requires that fluctuations around the av- 
erage energy contained in the jet above the cutoff and in 
the background medium are small. 

Let us stress at this point again the fact that while we 
can formally calculate the production of soft partons in 
a medium-modifled shower and flnd most of them still 
contained in a rather narrow cone (which would allow 
to get the total jet energy by making an angular cut 
without Pt cut and just subtract the uncorrelated back- 
ground), there is no reason to assume that soft par- 
tons or hadrons can be treated as part of a perturba- 
tive shower in the presence of a medium. On the con- 
trary, good evidence has been found in two particle cor- 
relations [33, UHl and three particle correlations [1^ at 
RHIC energies that the distribution of soft hadrons cor- 
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FIG. 6: (Color online) Differential jet shapes ^p{r) as a function of angle r given R = 0.7 for a 20 GeV shower initiating quark 
in vacuum and for three different in-medium scenarios (see text) for all final state particles (left panel) and a cutoff of Pt > 2 
GeV imposed (right planel). 
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FIG. 7: (Color online) Differential jet shapes ^(r) as a function of angle r given R = 0.7 for a 100 GeV shower initiating quark 
in vacuum and for three different in-medium scenarios (see text) for all final state particles (left panel) and a cutoff of Pt > 4 
GeV imposed (right planel). 



related with a hard trigger is very different in vacuum and 
medium. In essence, energy and momentum flow through 
soft hadrons is observed at large angles. This has been 
variously interpreted as jet-induced shock waves in the 
medium [s^, IH, [H, H^l . Given these finding, it is rather 
likely that a medium-modified jet is only a useful concept 
in terms of an energy flux in a relatively narrow angular 
region above a Pt cut, but that the soft dynamics is very 
different in vacuum and in medium. 

Thus, one needs to focus on the energy flux above a cut- 
off. If this is a quantity which does not fluctuate signifi- 
cantly, it may still be useful to estimate the total jet en- 
ergy. It is beyond the scope of this paper to discuss fluc- 
tuations of the uncorrelated background, but it is rather 
easy within a MC simulation to discuss fluctuations of 
the energy within a given set of cuts. 



A useful quantity in this context is the distribution of 
events dN/dE leading to a given energy E observed in- 
side the cuts, given that the jet was initiated by a parton 
with energy Eq. We show this quantity in Fig. [5] for a 
100 GeV quark as showe initiator in vacuum and for the 
three different in-medium scenarios for an angular cut of 
R < 0.6, the Pt cut of 4 GeV used previously and a 
more optimistic Pt cut of 2 GeV and the combination of 
the angular and the momentum cuts. Note that initial 
state effects as well as next to leading order pQCD dia- 
grams lead to a momentum imbalance between two back- 
to-back partons. This can effectively be treated with a 
Gaussian distribution of intrinsic fc^ which is added to 
the event, leading to the result that even in vacuum with- 
out any imposed cuts, the observed distribution of jet 
energy is never 6{E — Eg) but rather has a finite width. 
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FIG. 8: (Color online) Distribution of jets identified with energy E for a given set of transverse momentum or angular cuts for 
a 100 GeV quark as shower initiator. Shown are vacuum (upper left) and the three different in-medium scenarios RAD (upper 
right), DRAG (lower left) and FMED (lower right). 



It is evident from the figure that the medium leads to 
large fluctuations and that there is no good correlation 
between the energy detected inside the cuts and the ini- 
tial energy of the jet. The fluctuations are especially 
large in the radiative energy loss scenarios: Even if a 
rather optimistic cut of > 2 GeV is imposed, the en- 
ergy of the jet above the cut may be anything between 
25 and 80 GeV in the RAD scenario. The situation is 
somewhat better in the FIVIED scenario, but even then 
the fluctuations are much increased as compared to the 
vacuum case. 

These results would suggest that it is very hard to get 
a reliable estimate for the initial energy of a medium- 
modified jet without resorting to 7-jet correlations. On 
the other hand, given a set of cuts in R and Pt, there is a 
straightforward measure of the effect of energy flux out- 
side the cut (although this cannot quantify the full prob- 
ability distribution of finding energy outside the cuts) 
— the nuclear suppression factor Raa of jets (see also 
the discussion in [2a])- If the cross section for produc- 
ing a jet with energy E within a set of cuts R, Pt at 
rapidity y is {E; R, Pt) and in heavy-ion collisions 



^j^{E; R, Pt), then the nuclear suppression factor of 
jets is a straightforward generalization of the hadronic 

Raa as 



RfA'iE;R,PT) = ,J';t:(E"py (13) 



In a similar way as hadronic Raa < 1 indicates that 
one is biased to observe high Pt hadrons originating not 
from all events but only from those in which the medium 
effect was small can R^a < 1 be used to study how 
biased the sample of jets identified within given cuts is. 
However, hadronic Raa cannot simply be unfolded to 
find the probability distribution of energy loss 41], and 
for the same reasons jet Raa cannot in a straightforward 
manner be used to determine the probability distribution 
of energy deposition outside the cuts. 
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VII. CONCLUSIONS 

Identifying the precise nature of parton-medium inter- 
actions and hence also the relevant degrees of freedom 
in the medium is a major goal of the hard probes pro- 
gram within the physics of ultrarelativistic heavy-ion col- 
lisions. At present, different physics scenarios for this in- 
teraction, among them radiative energy loss, coUisional 
energy loss or a medium-induced drag force are all (more 
or less) compatible with the suppression pattern of high 
Pt hadrons as observed at RHIC (there are indications 
that the observed pathlength dependence in the suppres- 
sion of back-to-back correlations rules out a significant 
amount of elastic energy loss though [42.]). 
In this paper, it was shown that three such scenarios (two 
involving radiative energy loss, one a drag force) which 
lead to identical high Pt hadron suppression at RHIC 
[21] could be potentially distinguished using jet observ- 
ables such as the thrust distribution, the n-jet fraction 
or the jet shape. It should be stressed that even above a 
cutoff, the differences between the three scenarios appear 
substantial. 

The major obstacle to a measurement however would be 
the precise quantification of the initial jet energy. In 7-jet 
correlation measurements, this is rather straightfoward, 
however such measurements are difficult and have limited 
statistics. In contrast, a direct identification of medium- 



modified jets and a correction for the energy flow outside 
the cuts would appear difficult as the event-by-event fluc- 
tuations of energy outside the cuts are significantly larger 
in the medium for all scenarios under investigation than 
in vacuum. Note that the calculations presented here are 
also done for a fixed path. In a computation closer to the 
the experimental situation, the fluctuation of in-medium 
pathlength also has to be taken into account. 
It is tempting to try to devise a suitable analysis strat- 
egy for jet measurements in A-A collisions at LHC from 
these findings. This, however, would in all likelihood be 
premature: How a medium-modified jet at LHC would 
look like is at present rather model dependent, but an 
optimized jet finding strategy would need to make use of 
just this knowledge. It seems much more promising to de- 
velop strategies for jet finding and analysis iteratively by 
contrasting first data with qualitative theoretical expec- 
tations and then refine both the theoretical framework 
and the experimental analysis strategy as needed. 
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